// /opt/openmpi/bin/mpic++ -lm -o m mpiden.cpp
#include <mpi.h>
//    #undef SEEK_SET
//    #undef SEEK_END
//    #undef SEEK_CUR

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <math.h>

#include <iostream>
#include <fstream>
#include <iomanip>
#include <vector>
#include <queue>

using namespace std;
#include "general.h"
#include "mpiden_objects.h"
#include "mpiden_publicvars.h"
#include "mpiden_operations.h"
#include "mpiden_functions.h"

// masa por particula de la sim de 500Mpc es  7.1973e+10 masa solares
// masa por particula de la sim de 60Mpc  es  1.056e+8 masa solares

#define WRITERS 4      // Numero de files que se escribiran en el output
#define INTBOX   500   // 60
#define SNAPNUM  37    // 24
#define FILES    8     // numero de files por snapshot 
#define NPAQMAX 15000  // numero maximo de particulas por paquete de comunicacion

// Structuras a comunicar

int  totalsizeR,totalsizeL,maxsizeR,maxsizeL;
int main(int argc, char **argv)
{
  int i,j,ibox;	
  int iy, iz ,sendproc, recvproc;
  float numdens;
  float r0;
  FILE *pfout;
  long int kkk;
  char fileout[200];
  int blocklen[2]      = { 1, 1};
  struct ptcl tmptcl;
  
  MPI_Init (&argc, &argv);
  MPI_Comm_size (MPI_COMM_WORLD, &numprocs);
  MPI_Comm_rank (MPI_COMM_WORLD, &proc);

  // Crea la estructura Dosint
  interDer  = new struct face [2] ;
  delete[] interDer;


  r0 = 0.17;
  overhead = r0/2.;
  r0=0.17;
    
  nobj=92000;
  if(proc%2) for(i=0;i<100000;i++)
		for(j=1;j<10;j++)  r0=6*i/j;
  cout << r0 <<" soy proc " << proc <<endl;
  system("echo $HOSTNAME >> ttt.PASE"); 
  
  MPI_Finalize ();
  exit(0);
}

